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ABSTRACT 


Ocean  Acoustic  Tomography  requires  deep-ocean  moorings 
whose  horizontal  excursions  are  either  small  or  accurately 
measured.  The  present  study  rigorously  investigates  the 
former  case:  the  design  of  stiff  moorings  to  meet  any 
particular  horizontal  excursion  goal  (e.g.,  25  meters) 
under  two  typical  ocean  current-versus-depth  profiles. 
Moorings  are  considered  for  tomographic  transmitters  and 
receivers  at  depths  ranging  from  one  thousand  to  four 
thousand  meters.  Mooring  components  considered  include 
steel  sphere,  glass  ball,  and  syntactic  foam  buoyancy; 
jacketed  3  x  19  wire  and  electromagnetic  cable;  and  a 
realistic  (large)  battery  pack  for  the  acoustic  trans¬ 
mitter.  Kevlar  mooring  line  was  considered  and  rejected. 

The  basic  tool  of  this  study  is  a  well-verified  com¬ 
puter  program  that  simulates  mooring  motion.  Many  runs  of 
this  program  have  yielded  enough  data  to  make  plots  showing 
mooring  cost  as  a  function  of  excursion,  depth,  and  mooring 
type.  It  has  been  found  that  cost  and  excursion  are  very 
sensitive  to  the  current  profile.  For  moderate  currents 
(15  cm/sec),  placing  instruments  at  lesser  depths  (1000- 
2000  m)  is  expensive  ($50,000  to  $100,000  for  wire  and 
buoyancy  alone) .  Deep  deployments  are  much  less  expensive 
(about  $15,000). 

Given  in  appendices  are  studies  of  the  computer  simula¬ 
tion  of  ocean  internal  wave  currents,  power  systems  for  the 
long-range  acoustic  transmitter,  and  a  hypothetical  steel 
sphere  for  subsurface  buoyancy.  Also  given  are  tables  of 
mooring  component  characteristics. 
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INTRODUCTION 


Ocean  Acoustic  Tomograph} 


Ocean  Acoustic  Tomography  has  great  promise  as  a  means 
to  monitor  the  internal  "weather"  of  the  ocean.  By  measur¬ 
ing  acoustic  transmissions  between  many  moorings  over  large 
distances  (v  1000  km) ,  and  deconvoluting  by  computer, 
density  resolution  of  about  100  kilometers  should  be  possible. 
Furthermore,  depth  resolution  should  be  obtainable  by  dis¬ 
criminating  acoustic  multipaths  between  source  (Long  Range 
Transmitter,  LRT)  and  receiver.  For  a  detailed  discussion 
on  Ocean  Acoustic  Tomography  and  its  capabilities  see  Munk 
and  Wunsch  (1979) . 

There  are  two  obvious  sources  of  error  that  must  be 
strictly  limited  if  tomographic  acoustic  travel  times  are 
to  be  meaningful.  The  first  is  clock  drift'  between  source 
and  receiver,  which  must  be  Jess  than  about  25  milliseconds 
per  year.  The  second  is  motion  of  the  moored  sources  and 
receivers  -  any  change  in  separation  of  over  25  meters 
must  be  accurately  known  so  that  travel  times  can  be  cor¬ 
rected. 

This  position-keeping  requirement  can  be  met  in  two 
ways;  either  the  moorings  can  be  acoustically  tracked  by 
ocean  bottom-mounted  transducers,  or  else  the  mooring  can 
be  built  so  stiffly  that  its  excursion  will  rarely  exceed 
the  tolerable  limits.  The  latter  method  is  much  simpler 
than  the  former,  even  if  it  must  include  an  alarm  sensor  to 
warn  of  occasional  excessive  motion  caused  by  abnormally 
high  currents.  However,  because  of  unfavorable  scaling 
ratios,  stiff  moorings  can  be  massive  and  expensive  to 
build.  It  is  one  purpose  of  the  present  study  to  quantita¬ 
tively  say  how  massive  and  how  expensive. 
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B.  Background 


The  main  thrust  of  the  study  presented  in  this  report 
was  to  perform  preliminary  mooring  system  design  for  the 
first  large  area  experiment,  tentatively  planned  for  1985. 
Figure  1  is  a  functional  family  tree  illustrating  the 
moored  sensor  system's  role  in  relation  to  the  other  elements 
of  the  1985  ocean  acoustic  tomography  system. 

The  tool  used  to  determine  optimal  mooring  designs  is 
computer  simulation.  Over  the  past  few  years  the  authors 
have  created  a  set  of  computer  simulations  to  model  mooring 
mechanics  (CHHABRA  1973,  CHHABRA  1976,  CHHABRA  1977Aj. 

These  programs  have  been  verified  experimentally  (ChhABRA 
1974,  CHHABRA  1977^  CHHABRA  1979,  CHHABRA  et  al,  1974). 

For  the  present  study  the  simulation  is  used  for  two  things: 
(1)  To  determine  mooring  line  safety  factors,  not  only  in 
the  normal  final  deployed  state  but  also  during  an  anchor- 
last  deployment  free-fall.  (2)  To  predict  the  quasi-static 
3-dimensional  deflection  of  the  mooring  under  any  given 
current  vs.  depth  profile.  The  simulation  has  been  used  to 
design  various  source  and/or  receiver  moorings  that  meet 
the  tomographic  excursion  requirement  for  some  reasonable 
current  profiles. 

The  costs  of  these  moorings  have  been  tabulated,  based 
on  approximate  dollar  prices  as  of  late  1980.  Plots  have 
been  produced  of  cost  vs.  excursion  and  cost  vs.  source/ 
receiver  depth,  all  for  an  assortment  of  mooring  configura¬ 
tions.  The  details  of  this  work  follow  in  the  main  body  of 
this  report. 

As  a  by-product  of  this  study,  some  results  of  wider 
interest  have  been  obtained.  These  are  detailed  in  appen¬ 
dices,  and  include: 

(1)  how  to  synthesize  scientifically  reasonable 

current  profile  histories,  comprising  mean  uniform, 


1985  Ocean  Acoustic  Tomography  System 


FUNCTIONAL  FAMILY  1RLE 


mean  shear,  tidal,  inertial,  and  internal 
wave  components.  (Appendix  A) 

(2)  a  study  of  possible  power  systems  for  the 
Long  Range  acoustic  Transmitter.  (Appendix  B) 

(3)  a  detailed  set  of  tables  of  mooring  component 
characteristics.  (Appendix  C) 

(4)  a  theoretical  study  of  the  cost  and  weight 
of  ideal  steel  sphere  buoys  as  a  function 
of  working  depth.  (Appendix  D) 
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II.  MOORED  SENSOR  SYSTEM  GOALS  AND  SPECIFICATIONS 
A .  Scientific  Requirements 

For  tomographic  data  to  be  scientifically  useful, 
tomography  moorings  must  meet  exacting  requirements. 

A  tentative  definition  of  these  is  given  in  Table  1. 

The  most  challenging  of  the  specifications  are  the  two- 
year  lifetime,  the  real-time  data  telemetry,  and  the 
25-meter  horizontal  excursion  uncertainty  limits  on 
transmitters  and  receivers.  The  present  report  concen¬ 
trates  on  the  last  of  these  problems. 

The  lifetime  problem  is  addressed  by  the  explicit 
sizing  of  a  power  system  for  the  Long  Range  Transmitter 
(LRT) ,  and  by  implicitly  choosing  other  components  with 
a  two-year  lifetime  in  mind.  A  study  of  real-time  data 
telemetry  is  not  included  in  this  report;  a  hardware  de¬ 
sign  study  of  satellite  telemetry  via  a  pop-up-buoy, 
which  rests  safely  nested  atop  a  subsurface  moored  buoy 
except  during  brief  periods  every  few  days  when  it  sur¬ 
faces  to  transmit  data,  is  presented  in  DAHLEN  et  al. 

1981. 

An  obvious  approach  to  reducing  excursion  uncertain¬ 
ties  is  simply  to  stiffen  (increase  tension  in)  the  moor¬ 
ing  until  the  maximum  excursion  of  the  instrument  is  less 
than  the  allowed  uncertainty.  This  solution  is  generally 
possible  because  line  tensile  strength  increases  with  its 
diameter  squared  while  line  drag  increases  only  linearly 
with  diameter.  However,  because  line  weight  also  in¬ 
creases  with  diameter  squared,  a  law  of  diminishing  re¬ 
turns  eventually  takes  hold  to  forbid  arbitrarily  low  ex¬ 
cursions.  Of  course,  reliance  upon  this  stiff  mooring 
approach  presumes  a  priori  knowledge  of  at  least  the  gross 
features  of  typical  strong  current  profiles  at  the  mooring 
site.  It  may  also  be  necessary  to  provide  an  alarm  to  warn 
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1985  OCEAN  ACOUSTIC  TOMOGRAPHY 
SYSTEM  SPECIFICATIONS 


Source/Receiver  Horizontal  Excursion 


Maximum  Allowable 
Uncertainty 


Source 

Mooring 


±25  meters 


Source/Receiver  Vertical  Position 

Depth  Error 
Depth  Uncertainty 
Depth  Range 

Off-Bottom  Height  Range 


±25  meters 
±15  meters 
1000-4500  m 
>1000  meters 


Data  Telemetry 

Minimum  Interval  between  Messages 
Binary  Bits  per  Message 

Unattended  Lifetime  2  years 


Current  Profile 

Typical  1 
Typical  2 


see  Figure  2 
see  Figure  3 


TABLE  1. 


Receiver 

Mooring 

±25  meters 


±25  meters 
±15  meters 
1000-3500  m 
>1000  meters 


3  days 
3000 

2  years 
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of  the  occasional  occurrences  of  extreme  profiles  which 
push  the  mooring  beyond  the  tolerable  excursion.  A  brief 
survey  of  sensors  suitable  for  this  warning  is  presented 
in  Section  D  of  this  chapter. 

Another  approach  to  reducing  excursion  uncertainties 
is  active  position  tracking.  A  relatively  soft  mooring 
could  be  equipped  with  a  comprehensive  excursion  monitor¬ 
ing  system,  designed  to  measure  both  the  magnitude  and 
direction  of  excursions  to  better  than  25  meters  accuracy. 
The  most  likely  system  to  do  this  would  involve  acoustic 
ranging  between  the  source  or  receiver  and  transponsders 
on  the  ocean  floor  nearby.  Ideally  a  correction  for  each 
source-receiver  pair  would  be  computed  in  real-time  by  the 
receiver;  to  do  this  it  would  have  to  know  the  source's 
position,  sent  either  by  telemetry  or  else  coded  into  the 
long-range  acoustic  signal  itself.  Less  rigorous  correc¬ 
tion,  perhaps  using  average  position  and/or  post-process¬ 
ing,  etc.,  may  also  be  sufficient. 

B.  Current  Profiles 

The  mooring  motion  study  uses  two  profiles,  illus¬ 
trated  in  Figures  2  and  3 . 

Note  that  we  have  not  included  high  order  shear  or  any 
inertial  or  internal  wave  components,  since  a  crude  pro¬ 
file  serves  just  as  well  as  a  complex  one  for  comparing 
mooring  excursions.  If  we  wanted  to  simulate  time  his¬ 
tories  of  mooring  velocity,  then  a  full  complexity  pro¬ 
file  would  be  needed.  To  this  end  we  have  written  a  com¬ 
puter  program  to  synthesize  inertial  and  internal  wave 
currents,  using  a  rigorous  Garrett-Munk  spectrum  (see 
Appendix  A) . 

C.  Power  System  for  the  Long  Range  Transmitter 

This  item  earns  special  consideration  because  of  large 
mass  and  cost,  and  because  of  the  lack  of  engineering 
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attention  given  it  until  now.  In  contrast,  the  acoustic 
Receiver  and  Long  Range  Transmitter  are  under  development 
elsewhere  and  are  in  this  report  considered  merely  as 
black  boxes  of  given  size  and  weight. 

The  LRT  power  system  options  can  perhaps  best  be  sum¬ 
marized  in  a  family  tree,  drawn  in  Figure  4.  All  power 
system  options  assume  an  energy  requirement  of  85  kilo¬ 
watt-hours,  based  on  a  two-year  life.  Two  strawman  de¬ 
signs  for  a  transmitter  battery  are  sized  in  Appendix  B. 

D.  Alarm  Sensors 


In  the  stiff  mooring  approach,  each  mooring  would  be 
equipped  with  a  relatively  simple  excursion  monitor  to 
signal  an  alarm  during  the  occasional  high  current  event. 

On  a  receiver  mooring  this  alarm  would  be  telemetered 
along  with  the  data.  On  a  source-only  mooring  this  alarm 
could  be  sounded  by  the  LRT,  possibly  by  a  sequence  of 
omissions,  if  a  satellite  link  were  not  available.  Most 
likely  the  alarm  would  be  based  upon  a  3-day  excursion 
averaging  scheme.  Alarm  sensors  considered  are:  pressure, 
current  magnitude,  current  vector,  inclination  magnitude, 
inclination  vector,  and  tension. 

Pressure  Sensor 


A  simple  calculation  of  depth  excursion  for  a  given 
horizontal  excursion  and  mooring  line  length  gives 

Ah  =  x2/2l 

Where:  Ah  is  the  depth  excursion,  x  is  the  horizontal 

excursion  and  l  is  the  mooring  length.  For  a  2000  m  moor¬ 
ing  with  25  m  horizontal  excursion  Ah  is  16  cm.  This  seems 
impractical  to  measure  because  1)  the  resolution  of  the 
pressure  sensor  would  have  to  be  less  than  10  centimeters, 
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transmitter  (LRT) 


2)  averaging  would  be  required  to  take  out  the  tidal 
(%  1  m  amplitude)  fluctuations,  3)  the  procedure  requires 
knowledge  of  sensor  depth  at  zero  current  and  4)  drift 
errors  would  be  present  due  to  creep  in  the  mooring  line 
or  buoyancy  change  of  the  buoys. 

Current  Sensor 

A  simple  calculation  of  current  magnitude  (v)  required 
to  obtain  a  horizontal  excursion  (x)  gives: 

/  2Bx 

V  "  l 

where  B  is  the  buoyancy  (tension)  in  the  mooring  line  of 
length  t,  p  is  the  density  of  seawater,  CQ  is  the  drag 
coefficient  and  A  is  area  normal  to  the  current  direction. 
For  B  =  3000  lbs  =  13344  newtons,  x  =  25  m,  CQA  =  100  ft^ 
=  9.29  m2  ,  p  =  1025.9  kg/m^  ,  v  =  19  cm/sec. 

Magnitude  of  19  cm/sec  is  easily  measured  by  a  current 
meter.  An  unanswered  question  is  whether  magnitude 
measurement  at  just  one  location  gives  a  reliable  enough 
criterion  for  deciding  the  mooring's  real  horizontal  excur- 
excursion.  It  might  if  averaged  over  the  "right"  period, 
but  more  study  is  needed. 

Inclination  Sensor 

The  average  tilt  of  a  line  2000  m  long  and  having  a 
horizontal  excursion  of  25  m  is  given  as 

tan  6  =  25/2000  *8  =  0.0125  radians 

An  accelerometer  measuring  tilt  will  see  a  mean  of 
0.0125  g's,  which  is  a  good  size  signal.  The  dynamic 
inputs  to  the  accelerometer  would  have  to  be  either 
averaged  or  neglected  and  the  output  of  the  accelerometer 
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at  zero  inclination  would  have  to  be  known  with  the 
desired  accuracy.  The  most  serious  drawback  with  the 
tilt  sensor  is  that  it  could  not  be  placed  near  the 
top  of  the  mooring  at  the  LRT  or  receiver  where  the 
local  line  tilt  approaches  zero,  regardless  of  the 
average  tilt. 

Tension  Sensor 

Resolution  of  25  meters  is  not  practical  with  a 
tension  sensor.  For  a  2000  m  long  mooring  with  a  ten¬ 
sion  (buoyancy)  of  3000  lbs  the  change  in  tension  due 
to  a  horizontal  excursion  of  25  meters  is  0.23  lbs,  only 
1  part  in  13000. 

Only  the  current  sensor  seems  worthy  of  further  study. 
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E.  Alternative  Design  Approaches 

Various  approaches  were  studied  to  meet  the  performanc 
goals  of  the  tomography  moored  array  design.  For  each 
specific  design  approach  strawman  source  and  receiver  moor 
ing  designs  were  outlined.  Distributed  buoyancy,  discrete 
line  tapering  and  accommodation  of  the  full  range  of  de¬ 
sirable  source  and  receiver  depths  were  considered  as  sub¬ 
variations  within  each  strawman. 

The  first  approach  and  the  main  focus  of  our  design 
studies  is  outlined  below. 

Approach  A 

Separate  source  and  receiver  moorings. 

Stiff  source  and  receiver  moorings. 

A  separate  satellite  telemetry  mooring  placed 
alongside  each  receiver  mooring.  Intercommuni¬ 
cation  between  the  Receiver/Processor  (RP)  and 
the  Pop-Up-Buoy  (PUB)  is  accomplished  by  means 
of  a  two-way  acoustic  telemetry  link.  The  PUB 
is  a  satellite  communication  buoy  which  rests 
safely  nested  atop  a  subsurface  moored  buoy 
except  during  brief  periods  every  few  days 
when  it  surfaces  to  transmit  data.  It  is  des¬ 
cribed  in  DAHLEN,  et  al.  (1981). 

Ship-to-Long  Range  Transmitter  (LRT) ,  Ship-to- 
RP  and  Ship-to-PUB  two-way  acoustic  telemetry 
links  are  provided. 

Alarm  sensors  located  on  source  and  receiver 
moori ngs . 

PUB  can  transmit  but  not  receive  RF. 
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Source  Moorings (s) 


Alarm  can  be  sounded  only  by  omitting  part 
or  all  of  the  transmission  sequence. 

One  alarm  sensor  (most  probably  a  simple 
current  magnitude  sensor)  mounted  on  LRT . 

LRT  depth  1000  m  to  4000  m. 

Ocean  bottom  is  5000  m. 

High-power/low-drift  clock  can  only  be 
interrogated  via  the  Ship-to-LRT  acoustic 
link. 

Acoustic  Release  is  located  in  the  LRT. 

Conf igurat ions 

1.  Top  buoyancy  provided  by  a  syntactic  foam  sphere. 

LRT  and  the  top  buoyancy  at  the  same  depth. 
Lithium  primary  battery  placed  with  LRT. 

Tension  member  is  3  x  19  USS  wire  rope 
(jacketed) . 

2.  Same  as  1.  except: 

Top  buoyancy  provided  by  steel  sphere. 

3.  Same  as  2.  except: 

Steel  sphere  above  LRT  (at  much  less  pressure) . 

4.  Same  as  1.  except: 

Lithium  primary  battery  placed  upon  the  anchor. 
Electromechanical  cable  (also  the  tension 
member)  connects  the  battery  with  the  LRT. 

5.  Same  as  4.  except: 

Top  buoyancy  provided  by  steel  sphere. 

Same  as  5.  except: 

Steel  sphere  above  LRT  at  much  less  pressure. 


6. 


7. 


Same  as  1.  except: 

Top  buoyancy  provided  by  faired  glass  balls 

8.  Same  as  4.  except: 

Top  buoyancy  provided  by  faired  glass  balls 

Receiver  Mooring  (R) 

RP  depth  is  1000  m  -  3500  m. 

Ocean  depth  is  5000  m. 

One  alarm  sensor  (most  probably  a  simple 
current  magnitude  sensor)  mounted  in  RP . 

RP  contains  its  own  power  source. 

Dead-weight  anchor. 

Acoustic  Release  is  located  in  the  RP. 
Configurations 

1.  -  Top  buoyancy  provided  by  a  syntactic  foam 

sphere. 

3  x  19  USS  wire  rope  (jacketed)  used  as  the 
tension  member. 

RP  and  the  top  buoyancy  at  the  same  depth. 

2.  Same  as  1.  except: 

Top  buoyancy  provided  by  steel  sphere. 

3.  Same  as  2.  except: 

Steel  sphere  (top  buoyancy)  is  above  RP. 

4.  Same  as  1.  except: 

Top  buoyancy  provided  by  faired  glass  balls 

Satellite  Telemetry  Mooring  (T) 

Ocean  depth  is  5000  m. 

Top  subsurface  float  depth  is  45  m. 

Top  subsurface  float  is  a  steel  sphere. 


It 


USS  3  x  19  jacketed  wire  rope  is  used  as 
the  tension  member. 

Anchor  is  dead-weight. 

Acoustic  Release  is  located  under  the 
lowest  floatation. 

PUB  is  a  22"  ■*  sphere.  A  winch  is  located 
in  the  PUB  and  a  5/32"  <f  wire  (O.D.  3/16") 
is  used  as  the  tether  between  the  PUB  and 
the  top  subsurface  float. 


Configurations 


1.  All  buoyancy  (exclusive  of  PUB)  is  provided  at 
the  45  m  depth. 

2.  Buoyancy  is  provided  at  two  depths. 

One  at  45  m  and  the  other  deeper. 

Alternative  approaches  considered  but  not  studied  in 
detail  include: 


jroach  B 


Separate  source  and  receiver  moorings 
Stiff  source  moorings 

Combined  receiver/satellite  telemetry 
moorings.  Intercommunication  is  accomplished 
by  either  a  two-way  acoustic  telemetry  link 
or  hardwire. 

Excursion  tracking  of  receiver/satellite 
telemetry  mooring. 

Ship-to-LRT  and  ship-to-RP  two-way  acoustic 
telemetry  links  are  provided. 

Alarm  sensor  located  on  source  mooring. 
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Source  Moorings (s) 

Alarm  can  be  sounded  only  by  omitting  part 
or  all  of  the  transmission  sequence. 

One  alarm  sensor  (most  probably  a  simple 
current  magnitude  sensor)  mounted  on  LRT . 

LRT  depth  1000  m  to  4000  m. 

Ocean  bottom  is  5000  m. 

High-power/low-drif t  clock  can  only  be 
interrogated  via  the  Ship-to-LRT  acoustic 
link . 

Acoustic  Release  is  located  in  the  LRT. 
Receiver/Satellite  Telemetry  Moorings  (RT) 

Receiver  depth  1000  -  3500  m. 

Ocean  depth  is  5000  m. 

Receiver  contains  its  own  power  source. 

Dead  weight  anchor. 

Pop-Up-Buoy  is  nested  atop  subsurface  float 
at  45  m.  It  talks  directly  with  the  receiver 
module  (acoustic  or  hardwired)  or  to  a  ship. 

It  can  transmit  but  not  receive  RF. 

Subsurface  float  is  a  steel  sphere. 

Approach  C 

Combined  source/receiver  moorings. 

Stiff  source/receiver  moorings. 

A  separate  satellite  telemetry  mooring 
implanted  alongside  each  receiver  mooring. 
Intercommunication  is  accomplished  by 
means  of  a  two-way  acoustic  telemetry  link. 
Ship-to-LRT/RP  two-way  acoustic  telemetry 
link  is  provided. 

Alarm  sensor  located  on  source/receiver  mooring. 
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Source/Receiver  Mooring  ( SR) 


Top  depth  (Transmitter  depth)  1000  m  -  4000  m 
Receiver  depth  1000  m  -  3500  m 
Ocean  depth  is  5000  m. 

One  alarm  sepsor  (unspecified) 

/ 

Satellite  Telemetry  Mooring  (T) 

(Same  as  Approach  A) 


Approach  D 

Combined  source/receiver/satellite  telemetry 
moorings.  Intercommunication  is  accomplished 
by  either  a  two-way  acoustic  telemetry  link 
or  hardwired. 

Excursion  tracking  of  source/receiver/satellite 
telemetry  moorings. 

-  Ship-to-LRT/RP/PUB  acoustic  telemetry  link  is 
provided. 

Source/Receiver/Satellite  Telemetry  Moorings  (SRT) 

Source  depth  1000  m  -  4000  m. 

Receiver  depth  1000  m  -  3500  m. 

Ocean  depth  5000  m. 

Pop-UP-Buoy  is  nested  atop  subsurface  float 
at  45  m.  It  talks  directly  with  the  receiver 
module  (or  a  ship) .  It  can  transmit  but  not 
receive  RF. 

Subsurface  float  is  a  steel  sphere. 
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III. 


ANALYSIS 


This  section  presents  results  from  computer  simula¬ 
tions  of  mooring  excursions  for  single-point  moorings. 

The  cost  of  eacn  mooring  is  calculated  and  plotted  against 
excursion,  mooring  length  and  mooring  type.  Approach  A 
as  outlined  in  section  II  is  the  focus  of  this  analysis. 

The  cost  of  other  approaches  can  be  derived  from  the  cost 
of  approach  A. 

A.  Design  Assumptions 

1.  Wire  and  Rope  a)  U.S.  Steel  3  x  19  polyethylene 
jacketed  "Oceanographic  Wire  Rope"  of  constant  diameter. 
This  wire  rope  has  become  the  "standard"  for  oceanographic 
moorings.  Although  we  have  not  considered  alternatives  to 
USS  3  x  19  wire,  much  money  might  be  saved  by  switching  to 
a  more  competitively  marketed  "non-rotating"  wire.  Of 
course,  the  ocean  qualification  testing  such  wire  would 
require  is  rather  daunting.  b)  Industry-standard  poly¬ 
ethelene  jacketed,  single  conductor,  electromechanical 
armored  cable.  This  cable  is  not  precisely  torque-bal¬ 
anced  but  has  been  used  successfully  on  moorings,  c) 

Safety  factor  of  2  for  both  the  wire  rope  and  the  EM  cable. 
For  USS  3  x  19  wire  rope  used  for  subsurface  moorings  a 
safety  factor  of  2  has  proven  satisfactory  in  recent  de¬ 
ployments.  d)  Line  normal  drag  coefficient  of  1.4  (for 
flow  component  normal  to  the  line)  and  line  skin  friction 
drag  coefficient  of  0.01  (for  flow  component  parallel  to 
the  line)  . 

Kevlar  as  a  mooring  line  tension  member  was  ruled  out 
for  the  following  reason: 

For  a  given  working  load,  Kevlar  and  3  x  19  USS  wire 
rope  cost  about  the  same  but  Kevlar  has  a  larger  diameter. 
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Larger  diameter  means  larger  mooring  excursions.  On  the 
other  hand,  for  a  given  mooring  excursion,  Kevlar  mooring 
line  has  a  larger  working  load.  A  larger  working  load 
requires  more  buoyancy  and  even  larger  diameter  rope. 

These  make  the  Kevlar  mooring  significantly  more  expensive 
than  the  wire  mooring. 

In  comparing  Kevlar  with  wire  rope  different  safety 
factors  (breaking  strength  divided  by  working  load)  were 
assumed.  One  manufacturer  (Wall  Rope  Works)  of  Kevlar 
recommends  working  loads  for  non-critical  applications  to 
be  figured  at  20%  to  25%  of  average  breaking  strength. 

This  compares  unfavorably  with  the  50%  used  for  3  x  19  USS 
wire  ropes  in  recent  deployments.  For  a  given  large 
(greater  than  2000  lbs)  working  load,  Kevlar  (safety  fac¬ 
tor  -  4)  requires  a  larger  diameter  than  wire  rope  (safety 
factor  =  2) .  Larger  diameter  gives  proportionally  larger 
drag  and  larger  drag  means  much  larger  mooring  excursion, 
greater  than  proportionally. 

For  a  5000  lbs  working  load  a  suitable  Kevlar  rope 
(1/2",  average  breaking  strength  22,000  lbs)  costs  $1.21/ 
foot  (1982  prices) .  A  comparable  3  x  19  USS  wire  rope 
(5/16",  breaking  strength  10,300  lbs)  costs  $1. 02/foot 
(1982  prices).  But  because  the  wire  rope  weighs  0.097  lbs/ 
ft  more  in  water  (0.125  lbs/foot  versus  0.028  lbs/foot), 
it  requires  more  floatation  at  an  added  cost  of  $0. 43/foot 
(recent  price  of  $4. 40/lb)  to  obtain  the  same  average  ten¬ 
sion  in  the  wire  rope  mooring.  Hence,  the  effective  cost 
of  wire  rope  is  $1. 45/foot  compared  to  Kevlar's  $1. 21/foot. 
This  difference  was  insignificant  when  comparing  costs  of 
the  two  mooring  systems. 

Fairings  were  not  included  because  of  very  high  cost, 
deployment  difficulty,  and  sometimes  questionable  effec¬ 
tiveness  in  small  currents. 
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2.  Buoyancy  a)  Spherical  top  float.  The  float  is 
either  syntactic  foam  or  steel.  Steel  spheres  are  either 
models  marketed  by  ORE  or  are  our  own  hypothetical  designs 
(see  Appendix  D) .  Drag  coefficient  used  for  these  spheres 
is  0.05.  b)  Glass  balls  in  clusters  of  four.  Each 
cluster  is  packed  in  a  commercially  available  fairing 

2 

which  has  a  drag  coefficient  multiplied  by  area  of  3.2  ft  . 

In  general,  excursions  can  be  reduced  by  10-20%  if 
extra  buoyancy  is  allowed  to  be  distributed  along  the  wire, 
because  a  higher  average  tension  is  possible. 

3.  Instrumentation  a)  Receiver/Processor  package.  This 
is  an  8-inch  diameter,  6  foot  long  cylinder  and  weighs  300 
pounds  dry  and  150  pounds  wet.  It  includes  its  own  battery, 
b)  Long-Range  Transmitter.  The  package  has  three  cylin¬ 
ders:  1  of  1-foot  diameter,  20  feet  long,  flanked  by  2  of 
1-foot  diameter,  14  feet  long.  It  weighs  1500  lbs  dry  and 
500  lbs  wet.  Another  cylinder  houses  the  lithium  battery 
power  source  for  the  LRT.  Depending  upon  working  depth 
(1000  -  5000  m)  its  diameter  is  19-25  inches,  its  length  is 
45-47  inches,  its  weight  is  760  -  2620  lbs  dry  and  360  - 
1980  lbs  wet.  It  is  located  at  either  the  LRT  or  the  anchor. 
See  Appendices  B  and  C. 

Detailed  specifications  and  costs  of  the  mooring  lines, 
buoyancy,  instruments  and  anchors  used  in  this  study  are 
listed  in  the  tables  of  mooring  component  characteristics 
in  Appendix  C. 
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B.  Results 

On  the  following  pages  are  plots  of  cost  versus  excur¬ 
sion  and  cost  versus  depth,  for  various  mooring  types. 

These  are  given  independently  for  source  and  receiver  moor¬ 
ings,  and  for  the  two  current  profiles  described  in  section 
II.  B.  The  plots  represent  the  result  of  hundreds  of  runs 
of  computer  program  DYNOSUB. 

Of  particular  interest  are  the  plots  of  cost  versus 
depth,  with  the  horizontal  excursions  of  the  transmitter 
and  the  receiver  held  constant  at  25  meters.  These  plots 
provide,  in  summary  form,  a  direct  comparison  of  competing 
mooring  configurations.  A  full  description  of  each  configu¬ 
ration  is  given  above  in  section  II.  E.;  curves  on  the  plots 
below  are  keyed  to  these  by  a  code  signifying  Approach: 
Instrument:  Configuration.  For  example,  "AR3 "  denotes  a 
Receiver  mooring  using  approach  A,  configuration  3. 

Our  computer  runs  and  plots  explicitly  consider  only 
approach  A;  that  is,  separate  source  and  receiver  and  tel¬ 
emetry  moorings,  without  excursion  monitoring.  However, 
costs  of  moorings  using  other  approaches  can  be  inferred. 

For  instance,  the  cost  of  an  excursion-tracked  mooring  is 
simply  the  cost  of  tracking  equipment  plus  the  cost  of  a 
large-excursion,  soft  mooring  (given  on  the  cost  versus 
excursion  plots) .  The  cost  of  a  combined  source  and  re¬ 
ceiver  mooring  is  the  cost  of  a  source  mooring  plus  that 
of  the  buoyancy  needed  to  support  the  receiver/processor 
instrument  (because  the  drag  added  is  relatively  small) . 

C .  Conclusions  and  Observations 

1.  Most  of  the  drag  is  provided  by  the  line,  not  the 
buoys  or  instruments.  Therefore  the  thickness  of  the  line's 
jacket  should  be  minimized.  Excursion  will  be  very  sensi¬ 
tive  to  this,  and  relatively  insensitive  to  the  size  or 
weight  of  the  instruments  and  buoys. 
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2.  The  curves  of  cost  versus  excursion  are  almost  hyper¬ 
bolic;  cost  is  insensitive  to  excursion  when  excursions  are 
high,  but  excursion  is  very  insensitive  to  cost  when  excur¬ 
sions  are  low.  For  example,  for  profile  2,  all  excursions 
over  40-60  meters  are  cheap,  while  excursions  under  20-30 
meters  become  astronomically  expensive. 

3.  Buoyancy  materials:  glass  balls  are  the  cheapest  for 
deep  sources  and  receivers,  while  shallow  steel  spheres 
(about  450  m.  depth)  are  cheapest  for  less-deep  instruments. 

A  deeper  steel  sphere  (about  1000  m.  depth)  is  too  expensive. 
Syntactic  foam  is  also  more  expensive,  but  may  be  preferable 
to  glass  balls  for  intermediate  depths,  because  of  the  large 
number  of  balls  required.  For  example,  a  source  mooring  to 
have  27  meters  excursion  under  profile  2  at  2000  meters  uses 
240  17-in  glass  balls! 

4.  Electromagnetic  cable  between  the  LRT  and  its  power 
source  on  the  anchor  appears  to  be  the  most  economical 
approach  for  LRTs  deeper  than  about  1500  meters.  However, 
this  cost  advantage  diminishes  and  may  reverse  if  plain 
wire  is  used  (wire  that's  not  premium  priced  as  is  USS 

3  x  19)  . 

5.  The  bottom  line:  How  much  does  a  stiff  tomographic 
mooring  cost?  Assuming  25  meters  horizontal  excursion  under 
profile  2,  the  answer,  with  minimum  cost  configurations,  is: 


Instrument  Depth 
(meters) 

Receiver  mooring  without 
receiver/processor 
(1000  $) 

Source  mooring  with¬ 
out  LRT,  but  with 
battery*  (1000  $) 

1000 

108 

184 

1500 

95 

160 

2000 

60 

136 

2500 

39 

114 

3000 

25 

92 

3500 

16 

87 

4000 

10 

81 

*  includes  LRT  power  supply  for  $72,000. 

Note  that  the  receiver  moorings  cost  about  the  same  as  source 
moorings  minus  the  LRT  power  supply. 
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6.  Comparing  the  effects  of  profile  1  with  those  of  pro¬ 
file  2  shows  that  excursions  and  costs  are  very  sensitive 
to  the  current  profile.  Therefore,  for  any  deployment  site, 
it's  vital  to  have  a  good  prediction  of  the  particular 
currents  there,  so  that  a  suitably  stiff  mooring  can  be 
chosen.  For  currents  much  greater  than  profile  2,  moorings 
with  less  than  25  meters  horizontal  excursion  are  essen¬ 
tially  impractical  unless  the  instruments  are  deep. 
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Receive*'  Deo+h 

FIGURE  5 

Cost  comparison  of  mooring  configurations  for  the  acoustic  Receiver 
as  a  function  of  receiver  depth,  for  a  particular  current  profile 
and  horizontal  excursion. 
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Minimum  Receiver  mooring  cost,  as  a  function  of 
Receiver  depth  and  horizontal  excursion. 
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FIGURE  21 


Cost  vs.  excursion,  depth  for  mooring  design  AR2  under  current  profile  2. 
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APPENDIX  A 


Computer  Simulation  of  Ocean  Internal  Wave  Currents 


A  Fortran  language  computer  program  has  been  written 
to  create  synthetic  time  histories  of  North  and  East  (u  and  v) 
components  of  ocean  current.  This  is  accomplished  by  apply¬ 
ing  an  inverse  Fast  Fourier  Transform  (FFT)  to  spectra  defined 
as  follows: 

For  frequencies  greater  than  the  Earth's  inertial 
frequency,  the  amplitude  is  given  by  Equation  6-20  of 
Garrett  and  Munk  (1972) .  For  frequencies  less  than 
or  equal  to  the  intertial  frequency,  the  amplitude 
is  chosen  to  be  proportional  to  frequency,  providing 
a  decline  with  a  slope  of  two  on  a  log-log  scaled 
power  spectral  density  plot.  This  choice  is  motivated 
by  Figure  4  of  Garrett  and  Munk  reproduced  here  as 
Figure  1. 

Over  the  entire  frequency  range,  spectral  phase  is 
chosen  randomly,  independently  for  each  of  the  particular 
frequencies  in  the  discrete  spectrum  that  the  FFT  re¬ 
quires  as  input. 

The  spectra  for  North  and  East  components  are 
each  calculated  separately  in  this  fashion.  They  differ 
in  phase,  but  not  in  amplitude. 

After  the  spectrum  generation  and  inverse  FFT's  are  com¬ 
puted,  a  single  sinusoid  corresponding  to  the  tidal  current 
is  added  to  the  North  and  East  velocity  time  histories.  The 
tidal  amplitude  (uniform  over  depth)  is  assumed  to  be  5  cm/sec. 
(10  cm/sec.  peak-to-peak)  at  zero  phase  angle,  identically  for 
both  the  u  and  v  components.  Its  period  is  12.4  5  hours. 
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0. 


FIGURE  4  Bands  represent  the  observed  moored  and  floating  spectra  (l  eft,  see  Table  1 ) 
and  towed  spectra  (Right,  see  Table  2).  Dotted  bands  are  presumably  not  subject  to  fine- 
structure  noise.  Curves  are  the  derived  relations  (6.20,  21). 


(Taken  from  Garrett  and  Munk  1972) 


Figure  1 
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Equations 

1.  Spectrum  if  to  >  m  : 

Equation  6-20  of  Garrett  and  Munk  can  be  rewritten 

as : 


Ul 


1  + 


F  (<o)  =  8062 . 35(0^0)^  7 


tu^ 


V  1  -  -j 

CO 

2 

where  F  (to)  is  the  amplitude  squared, 
(meter^/hour^/CPH) 

(o  is  frequency  (CPH)  , 


1 

*  ~7 

to 


is  the  Earth's  inertial  frequency  (CPH), 


a>VB  is  the  Vaisala-Brunt  frequency  (CPH)  . 


Intuitively,  this  equation  can  be  divided  into 
three  parts: 

1.  A  semi-empirical  scale  factor, 

8062. 35to,,nujT  . 

Vo  X 

2.  An  ideal  degenerate  functional 
dependence, 

to 

This  produces  the  familiar  -2  spectral 
slope  on  a  log-log  PSD  plot. 

3.  A  slope-modifying  term, 
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This  approximates  unity  when  to  >>  to^.  , 
and  serves  to  produce  an  upsweep  in 
amplitude  as  u>  nears  to^  .  This  upsweep 
is  not  tame;  the  amplitude  is  infinite 
at  u  =  to  This  is  no  problem  so  long 

as  one  considers  the  continuous  spectra, 
because  the  area  under  the  peak  is 
finite.  However,  the  peak  could  cause 
a  large  phantom  excess  of  power  to 
appear  in  a  discretely  sampled  spectrum; 
If  a  sampled  value  of  to  happens  to  be 
very  close  to  ,  the  computed  power 

contribution  F^(oj)A(o  will  be  huge! 

As  insurance  against  this  hazard,  the 
computer  program  places  a  ceiling  of 
of  2000000  <o„_  on  the  value  of  F^(w). 

VD 


Spectrum  if  oo  ±  w  ; 

,  2 
(to)  =  2000000  -~ 


Earth's  Inertial  Frequency: 


(0 


I 


sin  (latitude) 
12 


4.  Vaisala-Brunt  Frequency: 

to  _  =  3  cycles/hour,  if  Z  >  mixing  layer  depth 
VB 

-z 

m  =  3  e  if  z  <  mixing  layer  depth 

VB 

where  Z  is  the  depth  in  meters. 
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Computer  Program  Notes 

The  present  program  is  an  extensively  revised  and 
corrected  version  of  one  written  by  Marvin  Gove  of  NORDA. 
(The  most  important  change  consisted  of  applying  the 

FFT  scale  factor  Aw  =  - - -  to  power 

total  time  of  data 

rather  than  amplitude . ) 

The  program  accepts  the  following  parameters  as  input: 

Title . 

Number  of  depths  (  separate  time  histories  to 
be  made  for  each) . 

Depths  (number  as  above) - 

Latitude. 

Time  increment  desired  for  u,  v  histories. 

Total  time  desired  for  u,  v  histories. 

Mixed  layer  depth. 

The  program  has  been  verified  by  "closing  the  loop", 
using  a  previously  written  Power  Spectral  Density  (PSD) 
program  to  re-create  the  original  spectrum  from  the  gen¬ 
erated  u  and  v  time  series.  At  each  stage  of  the  process 
the  total  power  has  been  calculated  and  observed  not  to 
change.  Figures  2  and  3  are  sample  log-log  PSD  plots  of 
North  velocity  histories  generated  by  the  program.  They 
conform  well  with  Figure  1,  the  ideal  spectrum  from  Garrett 
and  Munk.  Note  that  in  all  three  plots,  the  power  axis 
has  been  scaled  by  the  Vaisala-Brunt  frequency,  so  as  to 
make  the  curves  depth-independent.  They  are  still  latitude- 
dependent,  however. 

Figure  4  is  a  program  listing.  Figure  5  is  a  sample 
program  output. 
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LOG  PSD,  M  SQ.  /  HR  SQ.  /  CPH  SQ.  (NORM.  BY  V  -B  FREQ) 

-2.00  -1.00  0.00  1.00  2.00  3.00  4.00  5. 

n  i  j _ i _ i _ i _ i - 1 - 1 - 1 


T000 

jo  00 


LOG  PSD,  M  SQ.  /  HR  SQ.  /  CPH  SQ.  (NORM.  BY  V-B  FREQ) 

-2.00  -1.00  0.00  1.00  2.00  3.00  4.00  5.00 


U  FAKE  SPECTRUM  FROM  GMUNK A ,  100M  DEPTH. 
3  RUNS  OF  1024  SPACED  200.  28.69  DEG. 
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c 

c 

c 

c 

c 

c 

c 


THIS  r.-'O"-'."  CCMFU res  L-'SCTFELATED  KCRIZCMTAL 
DUE  TO  INTERN'  L  U'.'.E  MOTION  of  LVALUATIf  3  T;  E 
Fcv?:r»  tram  — m  cr  the  srrcii  j:, 

am  i>.rv*  ■_  .ti*i  ere  i  celcn  t.’E  mixed  live?  zs 
yaisala-:*;.  .r 


velocity  ccnfcnlmts 
in/eese  FAST 
;r  c-iveh  def;::. 

ASF'  -m  re?  THE 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LA.". LEI  L) 
UMAX 
XL  AT 
TI".X 
OUT 


Kz'.ZZP  IN  ICAA  FCMAT. 

MIX"'?  CF  CEFTH3  (10  tiAXIM'JM ) . 

i.atiti  ce  <c:smr.3>. 

mximum  Tins  iirecxesi. 

t:- increment  iseccncsi  re?  rrr. 


zmixl  =  mixed  laies  eerrii  c meters i. 
z1  ji  =  cefths  <;::tl~3). 


CU”'JT  re?  EACH  DEATH  13  <  U*.X  SETS  CF: 

Tie: 

HORIZONTAL  VELOCITY  CCMECMEMTS  U  AND  V  <  METIPS. 'CEOCX'D  I . 

(irr.f’E  kmax  is  the  latscST  re  ::r  cr  z  that  is  •••=  tmax/delt.  i 


C  THIS  Fre5~i*;  HAS  SWITTEM  BY  JIM  SCHOLTEN  CF  CSDL  DL'Er'E-  MM -JULY 
c  issc.  it  is  a:.'  extensively  revise:?  versecm  cr  a  am  by 

C  M'RVEN  CEVE  CF  MCTCA. 

C 


ism  ooce 

1  1  4  C.  C  J 

:?C4 

zz'i  ::co 

»  <**v»  «  *  *  ^ 

_  -  w  V  V  / 

ISM  0003 

ism  coo  a 

ISM  C010 
ISM  C011 
ISM  001E 
ISM  0013 


REAL'S  RN3M3,  KIDTH 

REAL*4  FJCKOO?),  UISIOOI,  V(E192I,  Fie.c 

1  lASLEUO.  moil  VSFHCIO, 

2  LI  I  10  It  LCIIO.SI,  L3i  10  >S  I ,  MTEMT ,  NTEMT 

INTEGER *•  4  XM.'.X,  MNUI 

EC’JIVALENCE  (UMAX,  r  .Tt  1 1 ) 

C;.7A  PI,  OTC'-V 
1  3.145936,  .0174553  / 

WIDTH  =  1.CC0 

FTC 71  =  0. 

F73T2  =  0. 


CCMDUTE  C-AF! 


7T-KUNK  CCEFF ICIENT  IECJA7ICM  6-CCI  «#«« 


C 

E  =  2.  *  PI  *  l.E-5 
E  IS  OIUENSICNLESS. 

N7ENT  5  3. 

NTEMT  IS  CYCLES  PER  HOUR ,  G/.RRETT-M'JNK  NCPMALIZATICM  FA 
MTEMT  =  0.1EEE-3 
MTEMT  IS  CYCLES  PER  METER 

CCEFF  =  C.  4  E  *  NTEMT  /  ( PI  *  MTEMT  **  Z I 


TCP  FOR  FFCC! 


;,«,;■***»**«»  INPUTS  «-.**»*«»»< 
C 


Figure  4 
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TON  0014 
I r-J  0015 
ISM  0016 
is:;  coi7 


READ  (5,5)  ( LABLEl L ) »  L=1 , 12 ) 

5  FORMAT  (10A4) 

READ  (5,*)  JtlAX,  XL  AT ,  DELT,  TMAX, 
READ  (5,*)  ( Z( J ) » J  =  1 » JMAX ) 


ZMIXL 


C 

ISM  C01G 

WRITE  (6,6)  ( LACLE( L  )  >  L  =  1 , 12 ) 

ISM  0019 

6  FC"“AT  (1H1,  1CX,  12A.4) 

ISM  0000 

'.."ITE  (6,13)  J:'*:.  XL  AT ,  SELT,  TMAX,  ZMIXL 

ISM  0001 

15  FC-MAT  (5X,  ' MO .  DEPOTS  =  ',  ,13,  5X,  ' LAT  (DEG)  ■ 

1  EX,  'CELT  (SI  ='.  c10.1,  SX,  'THAN  (S)  = 

2  SX,  *  ZMIXL  <tl)  =’,  F10.3,  EX,/) 

ISM  0002 

W"ITE  (6,32)  CCEFF 

ism  c::s 

32  FCRMAT  ('  GARRETT-MUMr",  CCEFFICIEMT  CM  SQ./M? 

V. 

C* *********  COMPUTE  IMIRTIAL  FREO-TMCY. 

ITS'  COC4 

CMI  -  SIMM  XLA7*DTC9  )  /  12. 

c 

CMI  IS  CYCLES.  KOLR. 

ISM  0005 

WRITE  (6,31)  CMI 

ISM  0006 

31  rCRMAT  (/•  INERTIAL  FREQUENCY  IS  (CFH):',F10 

ISM  0007 

7W2PI  =  PI  *  2. 

ISM  0000 

Ct'12  -  CMI  *  CMI 

ISM  COER 

V* 

KMAX  =  INK  TMAX  /  CELT  ) 

ISM  C030 

DO  24  I  =  8,14 

ISM  C 0 31 

J  =  2**1 

ISM  COSO 

IF  (J  .GT.  KMAX )  GO  TO  27 

ISM  CO' * 

24  CCMTINL'E 

ISM’  CO 35 

27  KMAX  =  J/2 

ISM  con 

SIMM  =  SRRT<  FLOAT!  KMAX)) 

ISM  0037 

KF  =  i:max  /  2  +  1 

ISM  COSO 

TT  =  DELT  *  KMAX  /  3600. 

c 

TT  IS  IN  i.C'JRS. 

ISM  0039 

CELCM  =  1  /  TT 

c 

OELCM  IS  CYCLES  PER  HOUR. 

C**»*».****i*  begin  depth  loop  *«****«».* 

ISM  0040 

CO  COO  J  =  I, JMAX 

c* 

>*•«««,<•  COMPUTE  LOCAL  VAISALA-DRUNT  FRECUENCY 

ISM  0041 

c 

IF  (Z(J)  .GT.  ZMIXL )  GO  TO  25 

ISM  0043 

VSFHl J )  =  3. 

:r:J  G0t4 

GO  TO  26 

ISM  C045 

25  VSFH(J)  =  3.  *  EXP(-Z(J>/1300.  ) 

c 

VSFH  IS  CYCLES/HOUR . 

ISM  0046 

26  CCMTINUE 

ISM  0047 

WRITE  (6,30)  Z(J),  V3FH( J ) ,  KMAX 

ISM  CC4S 

30  FC'MAT  ( ISO ,  I9X,  'DEPTH  (M)  =',  F10.3,  5X,  'VB 

1  F10. 3,  ’  KMAX  ,  110, 

2  //22X,  'TIME (HP 5) ' ,  3X,  ' U( HPS ) ' , 11X,  ' 

ISM  0049 

FU2( 1 )  =  0. 

ISM  0050 

c 

DM  =  0. 

*********** 


Figure  4  cont'd 
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C»»**«»*  LCAO  FU01KI  KITH  G-"~~FT-fVMK  SPECTRUM  rDD 

C  TO  FRO'.  IDE  FINITE  '.ALLIS  F C?  >  ROS'.;'.;  ISO 

IF  I  ;n 

IZ’".  7  H  i  N 

c 

THE  INERTIAL  r  RE  DJENO ! . 

ISM 

0C51 

DO  ICO  K  =  O.KF 

TC  { 

css; 

CM  =  CM  ♦  OELCM 

c:s3 

cm;  =  CM  *  CM 

ISM 

CC54 

IF  (CM  .GT.  C MI )  GO  TO  90 

ISM 

CCSs 

FLUFF  1  -  2CCCCOC.  «  0 M2  /  CMI2 

C 

fU2(K)  IS  METERS  SOL""E0  /  HOUR  SD'JAREQ  / 

CPH. 

*C‘* 

CC57 

GO  TO  95 

::.i 

ccr*$ 

90  C"R\T  =  CMI2  /  CM2 

ISM 

0:59 

FUOIKI  s  CGEFF  *  CMI  «  ( 1 . .CMRAT )  /  (  SCRT(  1 . -CNRAT  IT  M2  ) 

ISM 

CCoO 

IF  (FUOIKI  .GT.  2. Cot  FL'OC:)  =  2.E6 

ISM 

ccu; 

FU:(K 1  =  FU0<K  )  »  VBFHI J  ) 

IS  S' 

CC63 

«5  CONTINUE 

c 

WRITE  CF.M  AND  AMPLITUDE**:  FCR  DEEL'GGIND. 

c 

WRITE  (6,51)  CM,  FU2(K) 

c 

51  FORMAT  (  2M5.6  ) 

ISM 

00*4 

100  CONTINUE 

LOAD  U  AND  V  WITH  RANDOM  NUMBERS  BETWEEN 

0  AND  1,  TO 

DE  USED  IN  ASSIGNING  FHASES. 

ISN 

0065 

DO  50  K  =  l.KF 

ISM 

0066 

U(K)  =  O.SEO  +  RNDNR( WIDTH) 

ISM 

CC67 

V(K)  =  O.SEO  ♦  R.NDURf  WIDTH ) 

ISM 

0C63 

50  CONTINUE 

ASSIGN  RANDOM  PHASE  FACTOR  AND  GIVE  RCQ’JI 

RED 

C 

SYMMETRY  TO  FEAL  AND  INAGINASV  COMPONENTS 

CF 

c 

U  SPECTRUM. 

ISM 

0069 

Ir 

CO  i;o  K  =  l.MMAX 

ISM 

C070 

IF  (K  .GT.  KF 1  GO  TO  115 

ISM 

cot; 

PHASE  =  TWOPr  *  U(KI 

c 

THE  2  FACTOR  BELOW  IS  TO  CCMPc.MSATE  FCR  SPLITTING  THE  AMPLITUDE 

c 

AT  EACH  FREQUENCY  INTO  TWO  PIECES  OF  1/2  ANF 

LITL'CE, 

c 

A3  REQUIRED  BY  THE  FFT. 

ISN 

0073 

FU  =  SRRT  (  FU2!K)  /  2.  /  TT  ) 

ISM 

0074 

F(1,K)  =  CCS( PHASE )  *  FU 

ISM 

C075 

F(2,K)  =  SIN( PHASE )  *  FU 

IC'J 

0076 

GO  TO  1CI 

C077 

115  KI  :  2  «  KF  -  K 

ISki 

0376 

F(1,K>  =  F(  1  ,KI ) 

ISM 

C  0  •  9 

F(2,K)  =  -F(2,KI) 

c 

THE  ;  FACTOR  EELCW  IS  TO  COMPENSATE  FCR  SPLITTING  THE  AMPLITUDE 

c 

AT  EACH  FREQUENCY  INTO  TWO  PIECES  Or  1/2  AMP 

’LITUDE, 

c 

AS  REC'JIP.ED  EY  THE  FFT. 

IS.M 

0080 

121  PTOT1  -  PTOT1  *  0.5  *  (F(1.K)**2  +  F(2,K)**2) 

*  2.0 

ISM 

0081 

120  CONTINUE 

Lr 

r 

WRITE  (6,167)  ( (F(I,JJ), 1=1,2), JJ=1, 2043,3) 

c 

WRITE  (6,167) 

ISN 

0032 

c 

c 

167  FGRMAT  (3E15.6) 

(;*»**,,**  BRING  IN  INVERSE  EFT  TO  CONFUTE  U  FTCN  ITS  SPECTRUM. 


Figure  4  cont'd. 
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ISN  0033 


ISM  0034 
ISN  CCJ 5 
ISM  OC'o 
ISM  0037 
ISM  0C33 


ISM  CCS9 
ISM  CS50 
ISS4  CC72 
ISM  0073 
ISM  CC74 
ISM  03=3 
ISM  C3«o 
ISM  C097 
ISM  CC -3 
ISM  0099 
ISM  0100 

ISN  0101 


ISN  0100 
ISM  0103 
ISM  0104 


ISM  01C5 
ISM  OIOS 
ISN  0107 


ISM  0108 
ISM  0109 

ISM  0110 
ISM  0111 
ISM  0113 
ISM  0114 
ISM  0115 
ISM  0116 
ISM  0117 

ISM  0118 


C 

CALL  FCL'RT  (F,  KM,  1,  *1,  1,  0) 

C 

C  WRITE  (6,1671  ( ( F( I > JJ  > , I-I *2 ) > JJ-1 ,2048,8 ) 

C  WRITE  (6,167) 

C 

CO  140  K  =  l.KMAX 
U!  K  )  =  F(1,K) 

FTOTO  =  PTOTO  +  U(K)««2 
140  CONTINUE 

PTOTO  =  PTOTO  /  KM AX 
C 

C  DO  FFT  AGAIN  FOR  TEST. 

C  CALL  FCURT  (  F  ,MM,  1  ,-1,1,0) 

C  WRITE  (6.167)  (  (F(I,JJ), 1-1*2), JJ  =  1, 0048, 8) 
C  WRITE  (6,167) 

C 

£****,,„*  PEPEAT  FOP  V  CCMFCNGNT. 


DO  143 

K  =  l.KMAX 

IF  ( K 

. 3T .  KF)  GO  TO  ’ 

.35 

THASE 

=  TUORI  «  V ( K ) 

FU  =  S 

SRT(  FJ2(  K  )  /  2. 

/  TT) 

F  ( 1 ,  K  ) 

=  CCS( FHASE )  * 

FU 

F  (  2 ,  K  1 

=  SIHC  PHASE)  * 

FU 

GO  TO 

145 

135 

H  -  2 

*  KF  -  K 

F(l.K) 

=  F(l.KI) 

F  (  2 ,  K ) 

=  -  r  (  2  ,  K.  I ) 

145 

CONTINUE 

CALL  FCURT  (F,  KN,  1,  +1.  1.  0) 

C 

04*>.»4  A00  M2  EAPOTPOPIC  TIDAL  COMPONENT . 

C 

DO  150  K  =  l.KMAX 
VI K)  =  F(l.K) 

TM2  —  0.05  *  SIM( 1.405057E-4  *  CELT  *  (K-l)) 

C  TM2  IS  METEPS/SECC'.'O.  FERICD  IS  12.45  HCU~S. 

C  CONVERT  U  A':0  V  TO  N/SCC  FROM  M/HCUP : 

UtK)  =  U!K>  /  3--C0.  ♦  Tt(2 
VI K)  =  V ( K )  /  3600.  +  TM2 
150  CCMTINUE 
C 

C*»*4«*4*»*  OUTPUT  #»**«*•*»* 

c 

DO  200  K  =  l.KMAX 
T  =  FLOATiK)  *  OELT  /  36C0. 

C  T  IS  HCURS. 

WRITE  (10)  (T,  U(K),  V( K ) ) 

IF  (K  .ST.  45)  GO  TO  200 
160  WRITE  (6,165)  T,  U!K),  V<  K  > 

165  FCP.'XT  COX,  3t  F12.5,  5:.')) 

^00  CC v»  T I ; H 

WRITE  (6,166)  ( PTOT1 ,  PTOTO) 

166  FORMAT  <//•  U  POWERS  EEFC  E  AMO  AFTEP  INVERSE  FFT  ARE 

1  2E15.6,  '  M**2  / 

STOP 


Figure  4  cont'd 
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TEST  Pl'N  HUNGER  A10,  204S  FOIN7S 

ND.  DEPTHS  =  1  LAT  (DEG)  =  25.690  DELT  C  S  >  =  =  60.0  Tf'.A  •'  (S) 

G\RRETT-MUN.<  COEFFICIENT  (N  S3.  /HP  ./CYCLE )  =  6362.33  Zril/L  (M) 

INERTIAL  FREQUENCE  IS  (CP.H):  0.C4C006 


DEPTH  ( M )  =  ICO. 030  VE  FP2Q  (CrK)  - 


tin:; hrs) 

U(f»FS  ) 

vir.rs) 

0.06567 

-0.032-9 

-0 .05726 

O.E-3333 

-0.:: :i9 

-0.C37C0 

o.sccro 

-0.02.Ci 

-o.c.'::9 

1 .  C.6u7 

-0.cr.i6 

p  ^  i  / r  •' 

1 . 3-^33 

-0.CJ7J7 

-0.CI  "22 

1.6C0C0 

-0.0::.:: 

C.C0 737 

1 . £36d7 

-D. 0-4191 

o.cc*:: 

2.13333 

-0 . 0-32-* 

C.OIC^C 

2.*"Z CO 

-0.CC255 

G  .  0 4  6 

2 .  6  ci  67 

-0.02200 

0 .  c . :  0 

2.53333 

-0.02350 

c.072  7.3 

3.CCC23 

0.00515 

0.C9;75 

3 . 4o637 

0.C2611 

0 . 

3.73333 

0.07225 

o-o^roi 

•♦.cccoo 

0.07253 

0 . 0  72  72 

4 .  £  •_  t  o  7 

0. 03309 

0.07279 

4.53333 

0.0762  3 

0.07250 

4.22CC0 

C. 05914 

o.o>;.c-t 

5.C6.-67 

0.02E56 

c.c.:-.? 

5 . 3^333 

”  0 . 0  C  3  3  6 

G  .  C  5 

5.5-0220 

-0.0121c 

0.025 40 

5. 2-:  667 

-0 .04973 

0.0C1 06 

6 . 13333 

-0. 03622 

-0.C3J3  2 

6.vc:co 

-0.0  .326 

-0.03220 

o  .  6 6 7 

-0.C3  “24 

-  c .  c .  r  j  3 

6 . 93333 

-0.C3237 

-0 . 242-.  1 

7.202:0 

-0.000-0 

-  0  .  C  3  4  6  3 

7. ‘-'>6  6  7 

0.01123 

-0.0-2.37 

7.73333 

C.C1673 

-0.02571 

S.CCCC0 

0.070=2 

-0.02235 

S. 26657 

0.02215 

-0.C-.273 

E  .53333 

0.02613 

-0.222=3 

s.socco 

0.C21T., 

-0 . 04222 

9.C6-67 

-0.CC201 

-0.01512 

9.33333 

-0.01=17 

-0.232-72- 

?.c  OCOO 

-0.023SS 

-0.CE-.4i 

9 . 1 1>6  37 

-0 .22115 

-0.CCE2-0 

10.13333 

-0.01209 

-0.C60I4 

10.42020 

0 . Cl  234 

-  C  .  C r  1 0 

10.66557 

0.02325 

-0.02  2-24 

10.03333 

C.C 4249 

0. 00-31 S 

11.2CC00 

0 .07170 

-0.00230 

11.45667 

0.10077 

0. 00  =  73 

11.73333 

0.10531 

-o.oc »:o 

lC.OOCCO 

0.1072-* 

O.O:  i40 

U  FC’-’CRS  BETOPE  AND  AFTER  INVERSE  FFT  ARE  :  0.7691G1E*C5  0.72-91  -22*05  K>*2  /  H‘ 


Figure  5 
A-13 


'0  *-«  •*>  cr  o  .►  r*-  cl  cv  -{>  .  t  o  o  to  f-.  ,-i  n  to  til  cj  ei  •'!  rfj  o  -j-  <*i  •?  cf> 


APPENDIX  B 


Power  Systems  for  the  Long-Range  Acoustic  Transmitter 


STRAWMAN  #1 

Lithium  Primary  Battery 
Lithium  Thionyl  Chloride  (LiSOCl^) 


Assumptions 


(a)  Energy  required  is  85  KWH  derived  as  follows: 

(1)  2-year  life,  200  watts  peak  electrical  power 
required  during  acoustic  transmission,  0.5  watts 
standby  power  at  all  times,  3  minute  acoustic  trans¬ 
mission  every  hour  over  a  24-hour  period,  then  a 
48-hour  rest.  All  of  these  parameters  may  be  eased 
based  on  future  sea  tests. 

(2)  A  precision  time  reference  (clock)  is  required. 
There  is  some  evidence  at  hand  that  a  suitable  one 
might  be  obtained  requiring  1  watt  or  less  (con¬ 
tinuously)  .  Assume  lwatt. 

Therefore:  AVG  PWR  =  (i)  (~)  (200)  +  1.5  =  4.84  WATT 

TOT  ENERGY  =  (4 ' 84  \q~6o  ^ — — —  =85  KWH 

(b)  Use  ALTUS  battery  stack  of  17-inch  diameter  cells 
(cells  are  available  in  various  capacities  and  thick¬ 
nesses  between  0.5  and  7.35  inches).  Average  energy 
density  is  180  WH/lb.  and  15  WH/inch^. 

(c)  Battery  contained  in  a  cylindrical  pressure  vessel 
formed  from  ASTM  A285  Grade  C  pressure  vessel  quality 


B-l 


hot- rolled  steel  plate.  Working  stress  is  27,000  psi. 
Assume  the  vessel  has  hemispherical  end  bells.  Vessel 
walls  will  be  so  thick  that  buckling  is  not  a  problem 
(see  analysis  in  Appendix  D) .  Cylinder  inside  diameter 
is  18  inches. 

(d)  Pressure  vessel  cost  is  approximately  3  dollars  per 
pound  (based  on  current  prices  and  quantity  production) . 

(e)  Battery  cost  per  ALTUS  estimate  on  5  June  1980: 

one  unit  :  $100,000 

10  units:  70,000  ea 

25  units:  40,000  ea 

Calculations 

Battery  volume  =  5700  cubic  inches; length  =  25  inches 
Cylinder  length  is  therefore  26  inches 
Battery  weight  =  470  pounds 

An  approximate  stress  analysis  gives: 


where : 

R  .  =  outside  radius 

c/s 

rcys  =  inside  radius  (9-inches) 

0  =  working  compressive  stress  (27,000  psi) 

P  H  seawater  pressure 

P  (psi)  =.  1.46  x  depth  in  meters 
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Results 


Item/Depth  (m) 
Pressure  (psi) 

Rc  (in) 

Rg  (in) 

Vessel  Wgt  (lbs) 

Total  Wgt  (lbs) 

Seawater  Displ 
(lbs) 

Total  Wgt  in 
Water  (lbs) 

Vessel  Cost*  (§) 

Total  Cost*  (S) 

Added  Cost  of 
Floatation  ($) 


1,000 

2,000 

1,460 

2,920 

9.51 

10.09 

9.25 

9.53 

290 

640 

760 

1,110 

400 

440 

360 

670 

1,000 

2,000 

71,000 

72,000 

2,000 

4,000 

3,000 

4,000 

4,380 

5,840 

10.74 

11.48 

9.83 

10.17 

1,050 

1,540 

1,520 

2,010 

500 

560 

1,020 

1,450 

3,000 

5,000 

73,000 

75,000 

5,000 

7,000 

5,000 

7,300 

12.34 

10.54 

2,150 

2,620 

640 

1,980 

7,000 

77,000 

0 

(useful  anchor) 


*  assumes  10  -  25  units 


Note 

Published  data  by  Mr.  McCartney  of  NOSC,  San  Diego  indicate 
that  ALTUS  batteries  have  performed  well  in  ambient  pressures 
up  to  2000  psi.  There  is  some  hope,  therefore,  that  a  much 
lighter  containment  vessel  might  be  used  at  depths  in  the 
1000  -  2000  meter  range. 


STRAWMAN  #2 


Alkaline  Primary  Battery 

i 


Assumptions 

(a)  Energy  required  85  KWH,  as  for  Strawman  #1 

(b)  Use  E95  "D"  size  alkaline  cells,  VSTART  @  1.5V, 

Vend  @  • 9V  providing  9.28  AH  (from  pg.  310  of  Ever- 

ready  Handbook)  @  70°F.  Derate  50%  for  32°F  service, 
thus  assuming  4.64  AH  @  1.2V  avg  =  5.57  WH/cell.  Each 
cell  is  4.5  oz  and  3.3  in3.  Therefore  cell  energy 
density  =  19.8  WH/lb  and  1.69  WH/in3. 

(c)  Pressure  vessel  similar  to  that  used  for  Strawman  #1. 

(d)  "D"  cell  cost:  $.65  ea.  Add  $.50  ea.  for  battery 
assembly. 

Calculations 


Total  battery  requires  15,300  cells  weighing  a  total  of 
4,290  lb.  and  requiring  50,400  in3  (excluding  wiring,  pack¬ 
aging  and  containment  vessel) .  This  assembly  seems  totally 
impractical,  but  we  shall  press  on,  nevertheless.  Considering 
the  large  number  of  discrete  cells  to  be  packaged,  assume  a 
packaging  efficiency  of  75%  and  thus  a  total  volume  to  be 
enclosed  of  67,200  in3  .  Assume  assembled  pack  weight  =  4400  lbs. 

Total  battery  cost  is  15,300  x  $1.15  =  $17,600.  There  is  no 
evidence  that  these  cells  can  function  at  ambient  seawater 
pressure,  so  the  container  must  be  a  pressure  vessel. 

Using  the  same  terminology  defined  for  Strawman  #1,  we  find 
the  volume  of  steel  in  the  pressure  vessel  is: 
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(cylinder  volume)  =  V. 


1 


-1 


V 


cyl 


V  .  (end  bells  volume)  =  |  i  r  ^ 

spn  3  s 


-1 


«-§> 


3/2 


where  Vbat  is  inside  cylindrical  volume  for  the  battery 


(67,200  in3) 


Apparently  the  smaller  r  the  smaller  the  weight.  For 

s 

example  we  find  (  L  =  cylinder  length) : 


1000 

meters 

5000 

meters 

rs(in) 

L  (in) 

Vcyl‘in3> 

vsph<in3) 

V  .  (in3) 
cyl 

VsPh(in 

9 

264 

7,900 

266 

59,000 

1,850 

12 

149 

7,900 

629 

59,000 

4,380 

15 

95 

7,900 

1,230 

59,000 

8,550 

18 

66 

7,900 

2,120 

59,000 

14,800 

Since  rs  has  only  a  small  influence  on  total  weight, 
choose  a  convenient  value  for  ease  of  manufacture  and  manage¬ 
able  length: 

r  =  r  =18  inches  ,  L  =  66  inches 
s  c 


We  find: 


Item/Depth  (m) 

1,000 

2,000 

3,000 

4,000 

5,000 

Pressure  (psi) 

1,460 

2,920 

4,380 

5,840 

7,300 

Rc  (in) 

19.03 

20.18 

21.49 

22.97 

24.67 

Rg  (in) 

18.51 

19.06 

19.67 

20.33 

21.07 

Vessel  Wgt  (lbs) 

2,810 

6,120 

10,070 

14,800 

20,700 

Total  Wgt  (lbs) 

7,210 

10,520 

14,470 

19,200 

25,100 

Seawater  Displ 
(lbs) 

3,760 

4,200 

4,720 

5,350 

6,120 

Total  Wgt.  in 
Water  (lbs) 

3,450 

6,320 

9,750 

13,850 

18,980 

Vessel  Cost*  ($) 

8,000 

18,000 

30,000 

44,000 

62,000 

Total  Cost*  ($) 

26,000 

36,000 

48,000 

62,000 

80,000 

Added  Cost  of 
Floatation  ($) 

17,000 

32,000 

49,000 

69,000 

0  (useful 

anchor) 

*  assumes  10  -  25  units 

Note: 

The  total  cost  of  the  alkaline  battery  system  is  comparable 
to  or  greater  than  the  cost  of  the  lithium  battery  system  except 
at  depths  less  than  2000  m. 
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APPENDIX  C 


Tables  of  Mooring  Component  Characteristics 


Details  about  components  used  for  tomography  moored 
systems  construction  are  listed  in  Tables  C.l  through  C.8. 
The  components  listed  are  classed  as  lines,  buoyancy, 
instruments  and  anchor.  Characteristics  given  are  type, 
size,  shape,  weights,  breaking  strengths,  drag,  working 
depth  limit,  cost,  and  source.  The  data  presented  has 
been  compiled  by  the  authors  during  a  period  from  late 
1979  to  early  1981,  mostly  from  manufacturers'  specifi¬ 
cations  . 
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Table  C. 3 


BUOYANCY 


DESCRIPTION  COST  $  SOURCES  REMARKS 


Table 


APPENDIX  D 


Design  of  Hypothetical  Steel  Sphere  for 
Subsurface  Buoyancy 


Consider  an  ideal  hollow,  thin-walled  sphere.  Three 
of  its  properties  are:  Weight  in  air  (W) ,  Gross  buoyancy  (B) 
and  maximum  working  depth  in  the  ocean  (Z) .  It  is  shown 
here  that 


W  a  B  Z  (D-l) 

For  the  case  of  steel,  this  equation,  and  the  conse¬ 
quent  one  for  net  buoyancy,  are  evaluated.  A  plot  is  made 
of  air  weight  and  cost,  versus  maximum  working  depth. 

Gross  buoyancy  is  given  by  the  product  of  volume  and 
seawater  density: 

B  -  |  -IT  r3  Pwg  ,  (D-2 ) 

where  r  =  sphere  outside  radius,  inches 

3 

p  g  =  seawater  density  =  0.03709  lbs/in.  (on 
w 

average) 

B  =  0.1554  r3  lbs.  (D-3) 

For  a  thin  sphere  of  thickness  t  , 

W  =  4  it  r3t  p  g  ,  (D-4 ) 

m 

3 

where  pg  =  material  density  =  0.28  lbs/in.  for  steel. 


by. 


The  maximum  working  pressure  of  this  sphere  is  given 


P  = 


2to 

r 


,  neglecting  buckling 
(see  below) 


( D—  5 ) 


where  p  =  pressure,  psi 

a  =  working  material  stress,  psi 

=  27,000  psi  for  a  practical  steel  (for  example, 
ASTM  A  285  Grade  C  Pressure  Vessel  Quality 


hot-rolled  plate) . 

D-l 


J 


Using  equation  D-5  to  substitute  for  t  in  equation 
D-4  gives, 

W  =  4  ir  r"  p 


.3  _  Pmg 


2a 


3  P  d 


2  a  p 


■  [f  * r3  v] 


Using  the  numbers  for  steel: 


m 


w 


|  =  0.0004194  p 

where  p  is  in  psi 

|  =  0.0006125  Z 

where  Z  is  in  meters. 


Net  buoyancy  is: 


BN  “  B  -  W 


=  W 


( _ 1 

^  0 . 0006 


0006125Z 


-  1 


So: 


W  _  Z  __ 

Bn  1633  -  Z 

This  is  a  useful  result 

Quoted  prices  of  ORE  steel  sphere  buoys  between  56 
and  22  inches  diameter,  vary  between  $2.50  and  $5.00  per 
pound  dry  weight,  respectively.  We  assume  buoys  would  be 
purchased  in  quantity  in  a  few  large  sizes. 


(D-6) 


(D-7) 


(D-8) 


(D-9) 


D-2 


Thus:  $/W  51  2 

or 

$/Bn  51  22/  Cl 633  -  Z) 

where  is  in  pounds 

The  last  two  equations  are  plotted  in  Figure  D-l. 


(D-10) 


Buckling 


The  buckling  pressure  of  a  thin-walled  sphere  is 
given  by 


P  = 


0.365  E  t 


(D-ll) 


For  steel  E  =  29  x  106  psi. 

From  equations (D-5)  and  (D-ll)  we  see  that  buckling  resistance 
governs  sphere  size  at  pressures  below  275  psi  (188  m)  ,  and 
allowable  working  stress  governs  sphere  size  at  higher  pressures. 
In  our  application,  depth  is  always  greater  than  188  m.  Lesser 
depths  are  not  considered  in  Figure  D-l. 
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maximum  WORKING*  DEPTH 
(meter 

Figure  D-l  Steel  Sphere  Buoyancy 
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